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Abstract 

The effects of fluctuating boundaries on a superposition state of a quantum particle in a box 
is studied. We consider a model in one space dimension in which the initial state is a coherent 
superposition of two energy eigenstates. The locations of the walls of the box are assumed to 
undergo small fluctuation with a Gaussian probability distribution. The spatial probability density 
of the particle contains an interference term, which is found to decay in time due to the boundary 
fluctuations. At late times, this term vanishes and the quantum coherence is lost. The system 
is now described by a density matrix rather than a pure quantum state. This model gives a 
simple illustration of how environment-induced decoherence can take place in quantum systems. 
It can also serve as an analog model for the effects of spacetime geometry fluctuations on quantum 
systems. 
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I. INTRODUCTION 



The theoretical study of quantum systems is usually performed by assuming that for all 
practical purposes they can be considered as isolated. However, realistic physical systems 
are never isolated. Instead, they are always immersed in an environment that continuously 
interacts with them jlj. As a result of such interactions, quantum superpositions tend to 
be suppressed and are not usually found at the macroscopic level. Exceptions occur in the 
phenomenon of superfluidity, and possibly in gravitational wave detection technology, where 
quantum effects can arise on macroscopic scales 0, 0]- The mechanism behind the loss 
of coherence between the components states in a quantum superposition has been called 
environment-induced decoherence (decoherence for short) [4H6|, and it has been accepted 
as the mechanism responsible to the emergence of the classical world from the quantum 
physics 0, 0]. Decoherence has been subjected to intense investigation in recent decades 
and plays important roles in several areas of physics, including quantum field theory, con- 
densed matter, quantum optics, quantum-chromo-dynamics and cosmology (for reviews, see 
Refs. (9)- 12 1 and references therein). 

Imposition of fluctuating boundary conditions on quantum fields in certain physical con- 
figurations was proposed in Ref . [l3[ . The motivation was to address the problem of infinite 
energy density appearing when fixed boundaries are assumed It was found that bound- 
ary fluctuations can lead to finite energy densities and a natural description of surface energy 
densities. 

In this paper, we study a simple model of a particle described by a superposition of two 
quantum states inside a potential well in one space dimension. Fluctuating boundaries are 
introduced as a mechanism modelling environment interactions. The averaged probability 
density is derived and its time evolution is investigated. It is shown that the net effect of 
the fluctuations is to suppress the interference terms between the component states of the 
particle wave function. As time goes on, the resulting probability density reduces to the 
case corresponding to a statistical mixture of states. The results indicate that this model 
can be understood as an example of decoherence. The effects of the fluctuating boundaries 
are similar to the effects of fluctuating or time dependent electromagnetic fields on quantum 
charged particles (lo| . 

It is well known that boundary conditions in flat spacetime can model certain aspects of 
the effects of spacetime geometry or gravity on quantum fields. For example, the vacuum 
energy of quantum fields in an Einstein universe can be derived as a Casimir energy, using the 
same techniques used to find Casimir energies in cavities in flat spacetime [l6j]. This suggests 
that fluctuating boundaries might be used as analog models for the effects of fluctuations 
of the gravitational field or of the background spacetime geometry. We will not pursue this 
idea in detail in the present paper, but leave it as a topic for future study. 

The simple quantum mechanical system of a particle in a superposition of states of a 
potential well is developed in the following section. In Sec. Ill, we model the interaction 
between the system and its environment by means of fluctuating boundaries, described 
by a normal probability distribution. The averaged probability density of the particle is 
derived and its time evolution is discussed. A numerical example is presented at the end 
of this section. Final remarks are outlined in Sec. IV, including a short discussion on the 
extension of the model to quantum systems in three space dimensions and to the case of 
electromagnetic waves in a cavity. 
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II. THE MODEL 



We begin with a non-relativistic particle of mass m confined in an infinite potential well 
whose width is denoted by a. For the sake of simplicity we restrict our analysis to the 
one- dimensional case. Suppose that the normalized state of the particle in an arbitrary time 
t is given by a superposition of the first two available states as 

lf)(x,t) = -j=[ll> 1 (x,t)+ll>2{x,t)]. (1) 

The eigenfunctions ip n (n = 1,2) are independent solutions of the Schrodinger equation 
under the conditions that these functions vanish on the boundaries at x = and x = a, 
yielding 

«M) = yf sin (^r) e ~^' ( 2 ) 

with u n = n 2 7T 2 h/2ma 2 . Each one of these eigenfunctions describes a stationary state, as 
the corresponding probability density \ip n \ 2 is time- independent. However, time evolution 
occurs when the particle is governed by a linear superposition of ip n (x, t), such as that given 
in Eq. ([1]). The probability density can be obtained as 

m^t)\ 2 = \ (|^| 2 + |^ 2 | 2 ) + W coswt, (3) 

where we define u = UJ2 — wi, which represents the corresponding Bohr angular frequency 
of the system. As we see, the time evolution of \ip(x, t)\ 2 is exclusively governed by the 
interference term between ipi and ip2- In fact, \ip(x,t)\ 2 oscillates between its maximum 
(|^i I + I^DV^ and minimum (|^i| — |?/> 2 |) 2 /2 values, as depicted in the down inset frame 
in Fig. [TJ for a particular numerical model. 



III. PROBABILITY DENSITY WITH FLUCTUATING BOUNDARIES 

Now we wish to investigate the behavior of this quantum system when interaction with 
the environment takes place. In order to model the interaction between the system and 
its environment, we allow the positions of the physical boundaries to fluctuate under the 
influence of an external noise. This can simply be implemented by allowing the width 
parameter a to undergo fluctuations around a mean value a [l3||. We set a = a(l+e), where 
the dimensionless parameter e is described by a Gaussian distribution as 

with 9 related to the width o of the distribution by means of o 2 = 1/26. The mean value of 
a arbitrary function G(e) over e is a linear operation defined by 

/oo 
G(e)f(e)de. (5) 
-00 

Notice that (e 2 ) — (e) 2 = cr 2 , the mean squared fluctuation of e. 
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There does not seem to be a meaning to averaging the wave function ip(x, t), as it is not 
directly observable. (See, however, Ref. 18j for a discussion of the possibility of measuring 
a wave function in the context of a weak measurements approach.) Here we are interested 
in studying averaged values of observable quantities. Particularly, the modulus squared of 
the particle wave function, Eq. ([3]), represents the probability density associated with the 
position the particle . The average over width fluctuations of this quantity can be calculated 
by using Eqs. flSJ) and yielding 



(6) 



Here the angular bracket refers to the average over positions of the boundaries, which was 
defined in Eq. (jSJ). In what follows, each term appearing in the above equation will be 
considered separately. 

The two first terms appearing in Eq. (JBJ) can be expressed, using Eqs. (EJ) and (JSJ), as 




-de 1 



0(1 



sin 



mix 



0(1 + e] 



de 



(7) 



where n = 1,2. We assume a narrow distribution, tr < 1, which means that only small 



values of e contribute in Eq. ((7j). First use the identity sin u — [1 — Re(e m )]/2. Next 
Taylor expand to first order in e inside the exponential, but only to zeroth order otherwise. 
The integrations may be performed using the identity 



e ±iZe-ee* d£ 



The result is 



) ~ — — — cos 

' a a 



2nirx\ 
a J 



exp 



-2n 2 Tt 2 a 2 x 2 
U 2 



(8) 



(9) 



Further, as the particle is confined (0 < x < a), and the small a approximation, irax/a 1 
is assumed, we have 

2 . o /7Z7TX 



sin 



(10) 



a \ a 

Thus, the probability density associated with the energy eigenstates does not change signifi- 
cantly when boundary fluctuations are introduced. In other words, if the particle is initially 
in an energy eigenstate state ip n (x, t), it will remain in this state and its probability density 
will not undergo appreciable time evolution when small boundary fluctuations are present. 

Next, we consider the last term in Eq. (El), which describes the interference effects occur- 
ring in the system. From Eqs. (j2J) and (jSJ), we obtain that 



(\MHcosut) = ^-{Ai + AU-Ai-AU 
+A 1 + A_ x — A 3 — A_ 3 ) , 



(11) 



where A^ is defined by 



At 



1 + e 



exp 



qmx iut „ 9 
± -; ~ - 9e 



o(l + e) (1 + e) 



(12) 
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and u = (jj 2 (a) — u)i(a) = 3h7i 2 /2ma 2 . Proceeding as before, assuming <x < 1 and using 
Eq. (|SJ), we obtain 



A; 



? exp 



iqn 
a 



x ± 



x exp 



-g 2 7r 2 
A6a? 



x ± 



aoot 
qn 



2acut\ 
qw J 



(13) 



After a finite time, t > 1/w, the x dependence in the real exponential becomes unimpor- 
tant, and A^ can be approximated by 



At ^ 



? exp 



x ± 



out 
qn 



-rt 2 



(14) 



where T = 2u a . Using this result in Eq. ( fill , we obtain 



(|^l||^>2| coswt) 



1 



COS 



(t)- 



COS 



3ttx 



x cos(a;t) e 



-rt 2 



(15) 



As one can see this term describes oscillations modulated by a factor which decays expo- 
nentially in squared time. The time scale for the onset of this decay is 



1 

t = - ■ 



(16) 



In the case of an electron in a potential well with a w lA, this time is of order t a ~ 10 17 s. 
However, once the decay begins, the characteristic decay time is 



1 



(17) 



which is longer by a factor of about l/er. 

Combining the results in Eqs. ( TTOl) and ( fT5l) with Eq. ([6]), we find that the average 
(\ip{x, t)\ 2 ) becomes 



i 

+ r 
a 



cos 



TlX 



COS 



37ra;\ 



cos(a;i:)e" 



-rt 2 



(18) 



As time passes, the last term in the above equation falls to zero for t ^> td- Thus, the net 
effect of the fluctuations is to kill the interference term. Neglecting the last term in Eq. (fi~8]) 
we obtain that 



<i^r> = ^ii 2 +hA 2 r) 



(19) 



which corresponds to a weighted sum of probabilities as occurs when a statistical mixture 
of states is considered. 
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The magnitude of the parameter a depends upon the source of the boundary fluctuations. 
However, it is useful to make a simple model in which the boundaries are massive particles in 
a harmonic potential well. The ground state of a quantum harmonic oscillator is described 
by a Gaussian with characteristic width 



Ax 



h 



2Mu Q 



(20) 



where M is the mass of the boundary, and uq is the angular frequency of the harmonic 
oscillator. Plausible choices for these parameters might be a nuclear mass and the vibrational 
frequency of a nucleus in a molecule or a crystal lattice: 



Ax « 0.01 A 



30 amu \ 2 
M ) 



10 



15 , 



U 



(21) 



With these choices, and assuming a ~ lA, we have a ~ 10~ 2 . 

The behavior described by Eq. (jTSJ) can be illustrated by a numerical example. For 
instance, let us study the time evolution of the averaged probability density when fluctuating 
boundaries are considered in a particular model with x/a = 0.7 and o = 0.01. In this case 
Eq. fl6]) can be integrated numerically. The result is depicted in Fig. [TJ Alternatively we could 
have used the approximate solution given by Eq. (fl8|) . which leads to an identical graph, 
confirming the approximation used in obtaining Eq. ( 1X8]) . As the figure shows clearly, the 




OJ t 



FIG. 1: (color online). The figure shows the behavior of the probability density associated with the 
state defined by Eq. ([T]) when fluctuating boundaries are assumed. As time goes on, interference 
effects between the component states of ip(x,t) are suppressed. For fixed boundaries, no suppression 
of interference is found, as shown in the down inset frame. We set x/a = 0.7 and a = 0.01. 

probability density (IV'I 2 ) oscillates around the mean value ((|^i| 2 ) + (\fa\ 2 ))/2, converging 
to this value when t 3> td- As anticipated, the net effect of the fluctuating boundaries of the 
potential well is to kill the interference effect between the two state components ipi and fa- 
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Finally, when no fluctuations are present, the usual stationary solution holds, as illustrated 
by the down inset frame in Fig. [TJ 

In order to have an estimate, consider again the case of an electron in a potential well 
with a ~ lA. As shown in Fig. [T]the oscillations in (\ip\ 2 ) are completely suppressed when 
uot = 200, that is, 10~ 14 seconds after the boundary fluctuations are turned on. 



IV. DISCUSSION 

We have studied a simple one dimensional model of a quantum particle confined by 
fluctuating boundaries. A coherent superposition of energy eigenstates is quickly converted 
into a statistical mixture, with the interference term being damped in time as e~ rt . This 
seems to be a simple example of decoherence in quantum systems. This behavior can be 
understood as a loss of phase coherence. A quantum particle in a box can be viewed as 
undergoing repeated reflections from the box walls. The uncertainty in position of the walls 
leads to uncertainty in the phase of the reflected wave, which accumulates in time. 

Although we considered only one space dimension for simplicity, the same features are to 
be expected for a quantum particle in a three dimensional box with fluctuating walls. The 
crucial feature is the appearance of an oscillating interference term, such as that in Eq. (j^D, 
whose frequency depends upon the position of the boundaries, which will also be the case in 
three dimensions. It is not necessary that all of the boundaries undergo fluctuations. Even 
if only one boundary is subject to fluctuations, while all others are fixed, this will introduce 
a dependence on the fluctuation parameter in the time dependence of the interference term. 
When averages over the fluctuations are taken, a decaying exponential in squared time will 
result, just as in Eq. ffTol) . 

The same comments apply to standing electromagnetic waves in a cavity with fluctuating 
boundaries. We could average the energy density of such waves over the position of the 
boundary, and obtain essentially the same results as for a quantum particle. If a single 
normal mode of the cavity is excited, then the energy density will not be significantly altered 
by small fluctuations, as is the case for energy eigenstates of a quantum particle. However, a 
coherent superposition of normal modes will have an interference term in its energy density 
which will decay in time in the presence of boundary fluctuations. Note that this effect has 
nothing to do with the damping in the electric field oscillations due to power losses in a 
cavity. Rather it reflects a loss of the definite phase relation between the superposed modes. 

In summary, the effects of boundaries with fluctuating positions is a simple model for 
decoherence in a quantum system. The effects of fluctuating boundaries can model external 
influences such as fluctuating electromagnetic or gravitational fields. 
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